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density modulation. The allowed solutions are saturated ferromagnetism (FM), a spiral 
spin density wave (SSDW) and a double spin density wave (DSDW). The DSDW phase 
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metallic except at multiples of a quarter filling. Near half filling the uniform SSDW 
phase is unstable against phase separation into a half-filled antiferromagnetic phase 
and a hole-rich SSDW phase. The dependence of the ground state wave number on 
chemical potential is conjectured to be a staircase. Comparison is made with higher 
dimensional Hubbard models and the J\ — J2 Heisenberg model. 
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1. Introduction 

The Hubbard Hamiltonian, originally proposed as a model of itinerant magnets, has 
lately gained new interest as a possible Hamiltonian for the cuprate superconductors. 
The Hamiltonian is 

H = H + Uj2^n ll (1) 

i 

where 

H = -Y.tiAsCjs- (2) 

ijs 

Here i and j are site indices, s =f , j is a spin index and U is the on-site Coulomb 
repulsion. We shall be considering the one-dimensional model with nearest-neighbour 
hopping, tij = t for = 1 and t — otherwise, and arbitrary band filling (0 < n < 2). 
One of the few exact results is the Bethe Ansatz ground state for this case (Lieb and 
Wu 1968). An approximate solution of the same model, especially one that breaks 
symmetries present in the exact solution, therefore requires some justification. 

The aim of the present work is to investigate spin structures in the Hubbard model. 
Analogous questions arise in frustrated Heisenberg antiferromagnets, as reviewed by 
Schulz et al (1994). The simplest example is the one-dimensional Hamiltonian 

H — J i Sj • Sj+i + J2 ^ Sj • Sj + 2 (3) 

i i 

with Ji, J2 positive (Majumdar and Ghosh 1969). The classical ground state is a 
spiral; in the limit J\ <C J 2 , the even and odd sites decouple and the ground state 
becomes degenerate. In the spin 1/2 case the peak in the structure factor moves from 
the Neel value 7r to 7r/2 with increasing J 2 , as in the classical case (Tonegawa and 
Harada 1987, Zeng and Parkinson 1995). Unlike in the classical case, dimerization (local 
singlet formation) occurs for J2/J1 of the order unity or larger (Haldane 1982). The 
dimerization is not visible in two-spin correlation functions, but is seen in singlet-singlet 
correlations of the form ((So • Si) (Si • Sj+i)). 

In certain situations the classical Heisenberg model possesses a line of degenerate 
ground states (in addition to the global rotational symmetry). The standard example 
in two dimensions is a square lattice antiferromagnet with nearest and next-nearest 
neighbour exchange interactions, with J x < 2J 2 . The classical ground state here 
comprises two interpenetrating antiferromagnetic sublattices, with no coupling between 
the Neel vectors. The three-dimensional face centred cubic antiferromagnet is similarly 
frustrated, with collinear spin density waves degenerate with a continuum of non- 
collinear double and triple spin density waves (see e.g. Long and Yeung 1986). 
These degeneracies could be lifted "by hand" by adding various terms to the classical 
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Hamiltonian, for example anisotropy or biquadratic interactions. Such terms are not 
needed; "ordering by disorder" stabilizes the collinear phase, the disorder being either 
quantum (Shender 1982) or thermal (Henley 1987) in nature. On the other hand, Long 
(1989) and Henley (1989) have shown that non-magnetic impurities favour islands of 
the non-collinear phase localized about the impurity. Both effects are due to the large 
transverse susceptibility of each antiferromagnetic sublattice to fluctuations in the other 
sublattice. We will see here how an itinerant model provides another mechanism for 
state selection. 

Similar spin structures arise in the Hubbard model; indeed, this reduces to the 
Heisenberg antiferromagnet in the half-filled two-sublattice large-C/ limit. We shall 
be investigating Hartree-Fock (hf) solutions of the one-dimensional model. The 
paramagnetic HF state is always unstable for U > in one dimension. Perfect nesting 
ensures that the band susceptibility diverges at q — 2kp. The system is therefore 
unstable towards a spin density wave at this wave vector for infinitesimal U. For 
finite U nonlinear effects can lead both to a shift in this wave vector and a distortion. 
We recognise that mean field theories in general, and HF calculations in particular, 
overestimate the tendency towards magnetic order (which the Mermin- Wagner theorem 
forbids in one dimension). However, the ordering may survive in the form of short-range 
correlations. 

Most recent HF studies of the Hubbard model relate to the two-dimensional case, 
with reference to copper-oxygen planes in cuprate superconductors. An exhaustive 
search of solutions of the HF equations is impractical, with global minima hard to 
obtain and in any case dependent on the boundary conditions. Some authors restrict 
consideration to collinear magnetization, following Machida and Fujita (1984). These 
found an exact solution for a one-dimensional model with linearized dispersion: a soliton 
lattice with a snoidal spin density wave. Most studies find coplanar spin textures in 
the two-dimensional Hubbard model; however, Chubukov and Musaelian (1995) find 
evidence for non-coplanar textures. Verges et al (1991) identify a large number of 
stable or metastable configurations. Spiral spin density waves (ssdw), with wave vectors 
varying continuously with U and n (Dzierzawa 1992), become unstable near half-filling. 
Here the HF ground state appears to be a collinear soliton lattice in which the holes 
are localized on walls between Neel-ordered domains (Fujita et al 1991, Ichimura et 
al 1992). The same conclusion follows from a fourth-order Landau expansion, valid 
for weak coupling, with coefficients determined from the electronic structure (Schulz 
1990). Nearest-neighbour Coulomb repulsion tends to stabilize the SSDW phase against 
this hole clustering (Hu et al 1994). One may identify various competing processes: 
firstly, a collinear spin density wave will present the electrons with a spatially varying 
potential, which may open a gap at the Fermi energy and stabilize the collinear wave. 
On the other hand, large U might favour the more uniform charge density in a metallic 
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SSDW. Competing ordered states, such as charge density waves and superconductivity, 
should also be taken into consideration (Bach et al 1994); however, for positive U in the 
absence of other interactions the SSDW is favoured (Eriksson et al 1995). Some of the 
above features can already be seen in the one-dimensional case discussed in this work. 

Early interest of course concentrated on the three-dimensional model (Penn 1966). 
There have been studies of spin structures in itinerant antiferromagnets, notably face- 
centred cubic 7-Mn, where single and multiple spin density waves are degenerate 
ground states of the classical Heisenberg model. Hirai and Jo (1985) show how fourth 
order terms in transfer integrals lift this degeneracy. However, spin density functional 
calculations reveal a very small energy difference between these structures (Crockford 
et al 1991). Here HF studies of the one-dimensional Hubbard model can illustrate the 
origin of the energy difference. 

Bach et al (1994) have studied a generalized unrestricted HF theory, in which 
the quadratic HF Hamiltonian allows particle-nonconserving terms and therefore treats 
magnetic and BCS states on an equal footing. They report a number of theorems 
concerning the ground state symmetry, although not much is known rigorously about 
the finite-?/ Hubbard model away from half filling. 

The next section covers the solution of the HF equations. The solutions are restricted 
to uniform states, in which the local density of states is uniform up to spin rotation. 
Such states fall in two two-parameter families: spiral spin density waves (ssdw) and 
double spin density waves (dsdw). Each is parameterized by a field amplitude and the 
nearest-neighbour angle. The band structure is calculated for both families, and the 
energy is minimized with respect to the parameters. Section ^ presents the resulting 
phase diagram: the DSDW is the ground state only in a narrow region near quarter 
(and three-quarter) filling. Near half filling, the homogeneous SSDW is unstable towards 
phase separation into a half-filled antiferromagnetic domain and a hole rich (or electron 
rich) domain. Finally, section f| discusses the physical significance of the dsdw and the 
stability of the states found. 

A brief report of some aspects of this work has recently appeared (Samson 1995). 

2. Computational method 

2.1. The Hartree-Fock approximation 

The unrestricted HF approximation minimizes ({Aj, Wi}\H\{A.i, Wi}), the expectation 
value of the Hubbard Hamiltonian ([!]) in the space of Slater determinants \{A.i,Wi}). 
These states are ground states of the non-interacting many-electron system in a spin- 
and site-dependent Hamiltonian 

iy H F({A 4 , w^) = h + J2(~ A i ■ s * + w i ■ n *)> ( 4 ) 
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with 

n i = c }.s C is (5) 

s 

^ = lY, C i a *t C it- (6) 

A st 

We take the number of sites N a —>■ oo and work in the canonical ensemble, with a fixed 
number n of electrons per site. The general problem is a minimization of a function with 
a large number of local minima and saddle-points, dependent on boundary conditions 
and differing little in energy, in a 4iV a -dimensional space. We shall therefore be far 
less ambitious and restrict consideration to uniform states, defined as those with a site- 
independent spin-projected local density of states, referred to the local spin quantization 
direction. Defining the tight-binding Green function in the usual way as 

[G l3 ] st (E) = (is\(E - H^r'tit) (7) 

where \is) is a one-electron Wannier orbital, we require the local tight-binding Green 
function to be of the form 

Gu(E) = a(E) + b(E) ei ■ a. (8) 

The only spatial dependence allowed is in the local magnetization directions, given by 
the unit vectors ej. All atoms are equivalent and there is no charge density modulation. 
The fields are then uniform in magnitude: 
1 

Wi = an d ^ = Aej. (9) 

The energy then becomes 

E HF = min VHe({Ai,Wi}) (10) 

{Ai,wi} 

where the minimization is subject to the constraints @ and 

VW({A 4 , Wl }) = ^ ({Aj, Wi} \H HF \ {A,,^}) + ^ + hjn 2 . (11) 

This is equivalent to the self-consistency condition 

Aj = 2U ({Aj,Wj} |Sj| {Ai,Wi}} (12) 
(so that in the HF solution the field is parallel to the magnetization). 

2.2. Spin density waves 

The uniformity condition (|8|) implies restrictions on the allowed directions, as we see 
from the expansion 

G U {E) = g H (E) - ^A 9iA E >j ' °9ji( E ) + ■■■ ( 13 ) 
1 j 
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where g%j{E) is the tight-binding Green function of the band Hamiltonian H (|2[). 
Consistency between equations (|T^|) and (|8|) for all energies requires the conditions 

(ei_ fc + e i+k ) || ei Vz, k. (14) 

The configurations are therefore coplanar, which implies vanishing torque on the local 
moments (Small and Heine 1984). (Non-coplanar configurations are allowed in two or 
more spatial directions.) This allows only two classes of configuration, as shown in 
figure 0: 

(i) Spiral spin density wave (ssdw) 

ej = (sin Qi, 0, cos Qi). (15) 

(ii) Double spin density wave (dsdw) 

e 2k = (-l) fe (0,0,l) 

e 2fe+ i = (-l) fc (sin0,O,cos0). (16) 

(We take the lattice parameter to be 1, and the spins to be in the xz plane.) The SSDW is 
specified by two parameters, the field amplitude A and the pitch angle Q, varying from 
(the ferromagnetic phase TTTT • • •) to tt (the antiferromagnetic phase Tit! • • •)• The 
dsdw is similarly specified by the amplitude A and an angle 9. It may be thought of 
as two interpenetrating antiferromagnetic sublattices with the staggered magnetization 
vectors canted at an angle 9, varying between (the collinear or dimerized configuration 
TTII • • •) and 7r/2, which coincides with the Q = n/2 SSDW T^l^ - • • •• I n the classical 
Heisenberg model the molecular field of one sublattice on the other vanishes in this 
phase, giving a two-dimensional manifold of degeneracies; any translationally invariant 
two-spin correlation function is independent of 9. The dsdw is also a superposition of 
two ssdw states at Q — ±tt/2. 



2.3. Band structures 

This restricted HF system is invariant under a subgroup of the symmetry group of the 
Hubbard Hamiltonian. Although translational and spin rotational symmetry are both 
broken, the Hamiltonian remains invariant under a simultaneous translation and spin 
rotation and Bloch's theorem still applies. We defer discussion of the stability of such 



configurations against further symmetry breaking to section [O. 

We calculate band structures of the HF Hamiltonian @) in the standard way. 
Following Korenman et al (1977), we transform the spin quantization axis from the 
z axis to the local axis = (sin 6^, 0, cos#j): 

c is = J2 U st(0i)d lt (17) 
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where 

U{9) 



cos 0/2 -sin 0/2 
sin 0/2 cos 6/2 



The transfer integrals couple the dj and dj bands, and the band structures reduce 
to the solution of second- and fourth-order secular equations for the SSDW and DSDW 
respectively. A little algebra leads to the two bands 



EssBw(k) = -2t cos(Q/2) cos k ± y/A 2 /A + At 2 sin 2 (Q/2) sin 2 k (19) 
(with Brillouin zone —it < k < tt) in the SSDW phase and the four bands 



E DS Dw(k) = ±yA 2 /4 + 2t 2 ± v /A 2 t 2 (l + sin^cos2fc) +4t 4 sin 2 2fc (20) 

(with Brillouin zone — tt/2 < k < 7r/2) in the DSDW phase. 

Figure § shows sample band structures, illustrating how dimerization opens gaps in 
the folded SSDW bands. The Fermi surface follows from the condition that the length of 
k space occupied is 2im. For both families the Fermi surface always consists of either 
0, 2 or 4 points. 

2.4- Energy computation 

The energy of these configurations, following equation flTT|) , is 

^hf(A, Q or 6) = -L £ ( / E(k)dk) + % + \un 2 . (21) 

2-k /-< \ Je k )<e v / 4(7 4 



bands 



For the SSDW energies ( |T9|) this is evaluated in terms of the elliptic integral of the second 
kind, E(<j),k) (Gradshteyn and Ryzhik 1980): 



/ E SS jyw(k)dk = —2tcos(Q/2)smk F 
J o 

p 2 sin fcp cos fcp 



A 
±— 

2 



where 



1+p 2 E (tan _1 [Jl + p 2 tan k F ],p/ J 1 + p 2 ) - 



1 + p 2 sin 2 fcp 



(22) 



p= (4t/A)sin(Q/2). (23) 

We evaluate the DSDW energies by numerical integration of the band structure ( p0| ) . 

For each point (n, U) in the phase diagram, we minimize the energies of the SSDW 
and DSDW states with respect to the parameters (A, Q) and (A, 6) respectively. Because 
of perfect nesting, a non-zero SSDW solution always exists in some Q interval for U > 0. 
The DSDW phase requires a little more care, as the energy gain is small and a non-trivial 
solution is absent in much of the phase diagram. 
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3. Phase diagram 

3.1. Energies 

Figure |3](a) shows the SSDW and dsdw HF energies for U = At and < n < 1. Because of 
particle-hole symmetry, the range 1 < n < 2 contains no further information and will not 
be discussed further; references to quarter filling also apply to three-quarter filling. For 
comparison, the top curve in the figure shows the energy of the paramagnetic solution 
(A = 0) and the bottom curve is Shiba's (1972) numerical solution of the Lieb-Wu 
integral equations for the exact ground state energy of the one-band one-dimensional 
Hubbard model. While the HF energies are not a good approximation to the exact 
ground state of this model, their relevance is discussed in section |4.2j . 

The SSDW ground state corresponds to saturated ferromagnetism for small n. The 
pitch Q then increases continuously from 0, reaching tt at half filling (n = 1). For small 
U, the nesting condition implies Q ~ 2/cp = vrra; as U increases, the SSDW moves to 
smaller Q, reducing double occupancy. The dsdw phase is stable with respect to the 
paramagnetic phase in a small range of fillings, and is only stable with respect to SSDW 
in a very narrow region near quarter filling (n = 1/2). The dsdw is most stable at 
exactly quarter filling, when the energy is minimized at the collinear dsdw 6 = 0. The 
angle 9 increases smoothly and monotonically (initially linearly) with deviation from 
quarter filling. 

The above discussion applies for the range < U < Ux ~ 5.58£; for U > U\ the dsdw 
state always has higher energy than the competing SSDW ground state (with Q < tt/2). 
Only for U ~ At can the dsdw region be seen in such a plot. It is difficult to find a 
dsdw solution numerically for small U. A fit in the range t < U < 3t suggests that the 
energy gain in forming a dsdw has the form exp(— U2/U), with U2 ~ 17t. 

Suzumura and Tanemura (1995) have recently reported results on the HF ground 
state and excitations of the quarter-filled Hubbard model. They consider a different 
class of configurations, namely collinear Q = tt/2 spin density waves. Such states are 
pinned to the lattice by a small commensurability energy dependent on the phase of the 
wave. Their ground state, with phase 7r/4, is precisely the dsdw (of energy — 0.681£ at 
U = At, as calculated here). At all other phase angles, a Q = tt charge density wave 
coexists with the spin density wave. Maximum energy (— 0.679t) occurs at zero phase, 

corresponding to a magnetic configuration f ■ | • This is still less than our minimum 

SSDW energy of — 0.674i. They find the value of U 2 /t to be Att^/2 w 17.8, compared with 
our numerical fit of 17. 



9 



3.2. Phase diagram 

Figure |3](b) shows the phase diagram in the (n,U) plane. For U > U c (n), indicated 
by the bold line, the ground state is saturated (strong) ferromagnetism (fm) with the 
upper band empty. The absence of unsaturated ferromagnetism is a consequence of the 
band-edge divergence of the one- dimensional density of states. The energy of this state 
(for n < 1) is independent of U, as there is no double occupancy: 

Ehf = tsmnn. (24) 

7T 

The phase boundary (which appears to be linear for small U), can be calculated 
by expanding the band structure (|T9| ) to 0(Q 2 ) and integrating; the line where 
<9 2 Vhf(A, Q)/dQ 2 vanishes is 

, . 27m — sin 27m , . 

UJn) = . 25 

n sin mi 

The limit n — > 1 gives U c — > oo. Thus, in the infinite-repulsion limit, FM is (as expected) 
stable for all n ^ 1. At half filling (n = 1) the stable state is AFM. Figure |] shows the 
energy of the uniform ground state (fm, SSDW or DSDW) as a function of n for various 
values of U . Where the uniform ground state is FM, the energy coincides with the 
U = oo curve (|24"D . DSDW is, as already mentioned, stable in only a small region of the 
phase diagram. All states are metallic apart from the AFM at half filling and the DSDW 
at quarter filling, where the Fermi level lies in a gap. 



3.3. Phase separation 

A curious feature of the SSDW energy is that it is not a convex function of n near n = 1. 
The uniform phase is therefore unstable to phase separation. Such an instability in the 
Hubbard model was proposed by Visscher (1974) and is analysed here by means of a 
Maxwell construction (Marder et al 1990, Arrigoni and Strinati 1991). The broken line 
E ps (n) in figure §(a) coincides with ^hf(^) at n = 1 and is tangent to the curve at 
n = ni(U), which is the boundary for phase separation. For n\ < n < 1 the uniform 
SSDW is unstable towards a state with volume fraction (n — n\) /(l — ni) of the n = 1 AFM 
phase and (1 — n)/(l — rii) of the n = n\ SSDW phase. The energy of the phase-separated 
state is 

E ps (n) = (^-^H F (l) + (l-^ HF (n 1 ) < 

1 — n\ 

The boundary is shown as a broken line superposed on the phase diagram in figure 0(b). 
To the right of the boundary the uniform phase is unstable towards phase separation 
between the AFM phase and the fm or SSDW phase at the boundary. For U > 8t 
(approximately), the only phases are FM and AFM. Andriotis et al (1993), in supercell 
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calculations on the one-dimensional Hubbard model with collinear moments, see a phase 
separation of just this form for large U. In one dimension the phase separation is an 
artefact of the HF approximation, as the exact energy is a convex function of n (Shiba 
1972). If phase separation does occur in the Hubbard model in higher dimensions, it is 
a consequence of purely short-range interaction and would be suppressed by long-range 
Coulomb repulsion. On the other hand, we do not know whether the domains are of 
macroscopic size. The instability might signal phase separation at a more local level 
- the migration of holes to antiferromagnetic domain walls in a soliton lattice. This 
possibility depends on the sign of the domain wall energy. 

Because of the kink in the energy at half filling there is no separation into n > 1 
and n < 1 phases. A similar phase separation does occur between the SSDW and dsdw 
phases, implying that the pure DSDW phase is only stable at exactly quarter filling. To 
avoid complicating the diagram, this is omitted from figure ^| but is discussed in section 
12. 



4. Discussion 

4-1. Interpretation of the DSDW 

The dsdw phase illustrates a number of physical phenomena. Firstly, it is an extreme 
case of a commensurate soliton lattice: an antiferromagnet with (uncharged) domain 
walls on alternate bonds. Secondly, the collinear DSDW is clearly stabilized near quarter 
filling by the opening of a gap at the Fermi energy (figure |2|). 
The HF energy can be expanded in the field: 

^hf(A, {e,}) = V HF (0, {e,}) + -^A 2 £ • e l+k + 0(A 4 ). (27) 

a ik 

The qualitative features of the phase diagram follow from a theorem concerning the 
dependence of susceptibilities on band filling (Heine and Samson 1980, 1983). Consider 
the moment expansion of the Green function 

Gij(z) = [(z - H)-% = b l3 z~ x + H l3 z- 2 + H lk H kj z~ 3 + . . . , (28) 

k 

where H is a tight-binding Hamiltonian. The energy difference between two 
configurations A and B changes sign r times as a function of band filling if the leading 
term in the difference of Green functions has the asymptotic form 

Tr£[4 A) W-Gi B) W]^- r ~ 3 (29) 

i 

for large z, where Tr is a trace over spin. A consequence is that the effective nearest- 
neighbour exchange interaction J± changes sign twice in the band < n < 2, while J 2 
changes sign four times. At quarter filling, J\ is small and ferromagnetic and J2 is larger 
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and antiferromagnetic. Longer-range interactions are smaller. Hence the dsdw region 
of the phase diagram of the Hubbard model corresponds to the dimerized region of the 
Majumdar- Ghosh model discussed in the introduction. 

The same theorem, applied to nonlinear susceptibilities, describes the state selection 
between collinear and noncollinear dsdw. The leading term in the ^-dependence of the 
Green function must involve a biquadratic term in the angles. Leaving out the terms 
irrelevant to the present argument and summing equivalent paths gives 

TrG e 00 (z) = . . . - -Tr (A ■ atoiAi ■ at 10 A ■ at QX A x ■ at w ) z' 9 + . . . 

= ...-lA 4 t 4 [(e -e 1 ) 2 -|e xe 1 |V 9 + --- 

= ... - ^AH A cos(26)z- g + .... (30) 

The energy difference between collinear and noncollinear dsdw states with equal A 
therefore has six zeros as a function of band filling (and a cos 2$ dependence). This 
agrees with the numerical results: the collinear phase is favoured near quarter, half and 
three-quarter filling, and the non-collinear phase at other values. Similar behaviour is 
indeed seen in calculations on 7-Mn (Long and Yeung 1986, 1987). The two-dimensional 
analogue seems to be the windmill configuration observed by Ichimura et al (1992) in a 
quarter filled square lattice. We note that the conclusion about the sign of the energy 
difference is not restricted to small A, where a Landau expansion of the form ( p7[) can 
be truncated at fourth order (Schulz 1990). 

4-2. Stability and relevance of solutions 

The solutions presented here (except at quarter and half filling) are metallic. In a one- 
dimensional system we would still expect an instability towards a non-uniform state 
that opens a gap at the Fermi energy. Auerbach and Larson (1991) find just such 
an instability towards a local increase in the spiral pitch in the t — J model. A similar 
instability is seen in the two-dimensional Hubbard model (Zhou and Schulz 1995). Phase 
separation is a further indication of this instability. If the filling is rational, the distortion 
would be commensurate; if it is irrational, an incommensurate distortion leading to a 
Cantor set spectrum is a possible outcome (Ostlund and Pandit 1984). Early studies of 
such spectra were more concerned with charge density waves, which, unlike the uniform 
SSDW, modulate the lattice potential. A similar modulation in the latter case requires 
a further distortion of the SSDW, leading to coexisting charge- density waves. This will 
result in a small kink in the E{n) curves as in figure |j. The true HF energy, according to 
the Maxwell construction, will be the convex hull of E{n). The ground-state Q, plotted 
as a function of chemical potential for constant U, would then be a staircase. The dsdw 
results and the small lock-in energies suggest that any energy gain will be small, and 
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the tongues in the phase diagram will be correspondingly narrow. However, discussion 
of the fractal properties of this phase diagram will take us rather far from the physics 
of the Hubbard model. 

The results reported here are consistent with known exact results for unrestricted 
generalized hf (Bach et al 1994). These authors discuss inter alia the HF ground state 
of the repulsive Hubbard model on a bipartite lattice. This is AFM for half filling, and 
FM in the U — > oo limit for all other fillings. They cannot obtain similar results in the 
interior of the phase diagram. However, the large gap and the work of Suzumura and 
Tanemura (1995) strongly suggest that the collinear DSDW is the true unrestricted HF 
ground state at quarter filling. 

It must be admitted again that HF calculations fail to give the correct ground 
state symmetry, and give poor ground state energies; quantum fluctuations destroy 
AFM and SSDW order in one dimension, but may retain short-range correlations of that 
form. The higher-order correlations that distinguish the DSDW would be harder to see. 
Agreement appears to be better in two dimensions where HF solutions are a useful leading 
approximation (Mehlig 1993, Mehlig and Fulde 1994). For two and more dimensions the 
space to explore, even for uniform configurations, is much larger and an exhaustive study 
would be more time-consuming. On the other hand, symmetry breaking is possible in the 
ground state, so that HF calculations may give a ground state of the correct symmetry. 
The paramagnetic phase will be stable for small U (except for the half-filled bipartite 
lattice), and the unrestricted HF energy E{n) will be a smoother function than in one 
dimension, as it will be more difficult to open a gap. The calculations here can then be 
considered as a toy model, nevertheless giving physically relevant predictions of spiral 
phases and a special phase at quarter filling. 

While the HF results give little information on the integrable one-band one- 
dimensional model, the approximation is more useful in the degenerate case (which 
is not integrable). Many-electron atoms may be modelled by an iV-band Hubbard 
Hamiltonian with inter-orbital Coulomb repulsion (U/AN)nf and Hund's rule term 
— (I/N)Sf, dependent on the total charge and spin on the atom. There is no reason that 
these should be governed by the same coupling constant for N > 1. In the case U ^> /, 
Coulomb repulsion will suppress the phase-separated and non-uniform configurations, 
thereby stabilizing the SSDW against phase separation. The HF results for uniform states 
are independent of N, and are the leading approximation in the limit N — > oo. 

The energies obtained here may also be useful for finite temperature properties. In 
a study of the thermodynamics, Samson (1989) fitted the SSDW energies to an extended 
spherical model, obtaining indications that the correlation functions fall more slowly 
with temperature than in the Heisenberg model. That approximation however cannot 
distinguish the DSDW from the Q = n/2 SSDW; the small energy gain suggests they would 
have little thermodynamic significance. 
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While HF solutions of the one-dimensional Hubbard model clearly cannot give 
definitive answers to questions about higher dimensions, they do provide a model in 
which many of the physical processes underlying state selection operate in a transparent 
way. 
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Figure captions 



Figure 1. Magnetization directions in (a) the SSDW and (b) the DSDW phases. 



Figure 2. Band structures for (a) the SSDW with A = t and Q = tt/2 (equivalent to 
the DSDW with A = t and 6 = tt/2) and (b) the DSDW with A = t and = 0. 



Figure 3. (a) HF results for the ID Hubbard model as a function of band filling n 
with U = At. The left-hand axis shows the pitch Q of the SSDW (bold dashed line). 
The right-hand axis shows the energies. From top to bottom these are paramagnetic 
phase (A = 0) (dash-dot line); SSDW phase (dotted); DSDW phase (full line); exact 
solution (bold line). The short dashed line on the right is the Maxwell construction, 
(b) The HF phase diagram for uniform phases (full lines). Between the dashed line 
and n = 1 the uniform phase is unstable to phase separation. 



Figure 4. Energy of the uniform HF ground state as function of band filling n for 
various values of U . The U = oo curve corresponds to FM. 
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